QED with a spherical mirror 
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We investigate the Quantum-Electro-Dynamic properties of an atomic electron close to the focus 
of a spherical mirror. We first show that the spontaneous emission and excited state level shift of the 
atom can be fully suppressed with mirror-atom distances of many wavelengths. A three-dimensional 
theory predicts that the spectral density of vacuum fluctuations can indeed vanish within a volume 
A 3 around the atom, with the use of a far distant mirror covering only half of the atomic emission 
solid angle. The modification of these QED atomic properties is also computed as a function of the 
mirror size and large effects are found for only moderate numerical apertures. We also evaluate the 
long distance ground state energy shift (Casimir-Polder shift) and find that it scales as (A/7?) 2 at 
the focus of a hemi-spherical mirror of radius R, as opposed to the well known (X/R) 4 scaling law 
for an atom at a distance R from an infinite plane mirror. Our results are relevant for investigations 
of QED effects, and also free space coupling to single atoms using high-numerical aperture lenses. 



PACS numbers: 42.25.-p, 12.20.-m, 37.30.+i 

Spontaneous emission and level shifts of atoms can be 
notably altered by placing them between mirrors. By 
modifying the electromagnetic mode structure interact- 
ing with the atomic electron [1] , one obtains a significant 
change in these quantum-electrodynamic (QED) atomic 
properties. Most experimental studies make use of high 
finesse cavities [2HZj to see the effects. Another way 
to change the properties of single emitters is to place 
other identical atoms close-by as originally propounded 
by Dicke |S]. To observe large QED effects in this case, 
the dipole emission patterns have to overlap, which re- 
quires the atoms to be very close to each other. Such 
effects were analyzed using two trapped ions [5J, but the 
Coulomb interaction between the ions restricted their dis- 
tance to a few microns. The interaction between two neu- 
tral atoms is however not overwhelmed by the Coulomb 
force. Using the large dipole moments of nearby Rydberg 
atoms localised in a dipole trap, entanglement between 
neutral atoms was recently demonstrated |10[ lllj . 

In general, an atom close to a single mirror already 
provides a very efficient way to investigate QED effects. 
The resonance fluorescence of a Doppler cooled Barium 
ion was reflected back onto itself in [12] , using a large nu- 
merical aperture (NA) lens and a mirror that was 30 cm 
away from the ion. In this experiment, the description of 
the interaction between the atom and the modified elec- 
tromagnetic field, or the mirror image, is very similar to 
the direct dipole-dipole coupling between two real atoms. 
Here, due to the high numerical aperture of the collection 
lens, the mode structure was altered significantly even if 
the mirror was many wavelengths away from the ion. A 
1% change in the decay rate was measured and found 
to be mostly limited by the collection solid angle and 
residual atomic motion. Such a system also leads to a 
vacuum-induced level shift in a laser-excited atom. This 
has been measured in |13j and found in good agreement 



with theoretical work [13] . 

A closely related field of research investigates the ab- 
sorption of photons from single atoms in free space. The- 
ory predicts that the best possible absorption between an 
incoming field and single atom arises when the incoming 
field matches the spatial atomic radiation mode [El [16] . 
Recent experiments have demonstrated substantial ex- 
tinctions from single molecules [T7lU9j . atoms [501 121] 
and quantum dots [22] in free space, thus showing the 
potential of free space coupling with high NA optics for 
fundamental investigation of light-matter interactions. 

The above mentioned studies make use of the inter- 
action of real photon with single atoms. There has 
also been a rapid increase in the number of experiments 
related to Casimir forces between dielectric materials, 
which are the result of the modification of the mode den- 
sity of virtual photons. Such studies are now being under- 
taken with an unprecedent level of precision (see for ex- 
ample [23] and references therein) . The comparison with 
the theory in this field is now reaching good agreement 
for some geometries, and over a wide range of materials. 
The possibility to use these measurements to gain a bet- 
ter control over nano-mechanical systems, and for precise 
tests of QED has been a major force driving this research. 
Although many geometries have been investigated theo- 
retically over the past decades [53], there are still many 
investigations relating to the sign of the force [25] , or the 
accuracy of the proximity force approximations [26], for 
estimating Casimir shifts of various materials. 

For atoms close to dielectrics, the modification of the 
ground state level shift (Lamb shift) yields the well 
known Casimir-Polder force [27], that was observed in 
[7J for a plane mirror geometry. The Casimir-Polder force 
was not reported nor calculated for single well localized 
atoms around complex opened 3D geometries. It is ex- 
pected that such investigations would also provide effi- 
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cient ways to test the behavior of vacuum fluctuations. 

Here, we demonstrate that a spherical mirror cover- 
ing half of an atomic dipole emission profile can fully 
suppress its spontaneous emission and excited state level 
shifts provided the mirror is close enough to allow tem- 
poral interference of the field amplitudes (Markov limit). 
We also calculate the shift of the ground state, the Lamb 
shift, as a function of mirror distance and found a scaling 
law that is more favorable than the plane mirror geom- 
etry for observing large shifts. Due to the development 
in the control of atomic motion |28) and mirror and lens 
fabrication [2"5H3T] , these effects are now within experi- 
mental reach. 

The modification of spontaneous emission and level 
shifts is first calculated using a one-dimensional model 
where the electron radiates along the mirror axis. We 
find full suppression of the vacuum fluctuations coupling 
to the atomic electron, even with a single mirror. The 
physical origin of this complete cancelation lies in the 
high spatial and temporal interference between the plane 
waves modes going to the mirror and the modes going to 
free space so that the density of vacuum fluctuations can 
reach zero around this idealized atom. Similar calcula- 
tions were performed by several authors to model more 
realistic scenarios (see for example [14, 32 and references 
therein) with the inclusion of free space vacuum modes 
that do not interfere with the mirror modes; therefore 
full extinction of spontaneous mission was not explicitly 
considered. We furthermore extend the one-dimensional 
calculations to a three-dimensional theory that goes be- 
yond the paraxial approximation, and demonstrate that 
the effective coupling to vacuum modes around an atomic 
electron can also reach zero within a volume A 3 around 
the focus of a spherical mirror. In the spherical basis, 
such an effect can be understood as the result of an inter- 
ference between even and odd spherical modes. Last, the 
long-distance ground-state energy shift (Casimir-Poldcr 
shift) is evaluated using the complete set of normal modes 
of the spherical mirror. We find that the Casimir-Poldcr 
shift scales as (X/R) 2 for a half mirror of radius R and 
atomic transition wavelength A, as opposed to the well 
known (A/i?) 4 scaling law for a plane mirror, where R is 
the mirror distance from the atom. 

There exists a wealth of studies on this topic. We 
would like to point out that, to the best of our knowledge, 
theoretical efforts concentrated on plane geometries [32T 
135] , on cavities in the paraxial approximation [HE], di- 
electric spheres [36H40] , or parabolic mirrors [41] . Closely 
related work was however done for a large spherical open- 
cavity by J.M. Daul and P. Grangier [32] • Strong en- 
hancement and inhibition of vacuum fluctuations were 
found with moderate cavity finesses and a full set of nor- 
mal modes were derived for the symmetric geometry. The 
case of a single mirror was found by setting the second 
mirror reflectivity to zero in a more general formula for 
the density of vacuum fluctuations in an asymmetric cav- 
ity. In section |III[ we present a different route towards 
finding the normal mode amplitudes in the single mir- 



ror geometry, that we further use for the Casimir-Polder 
shift calculation. 

We would like to stress that all the calculations are 
performed in the limit where the mirror is many wave- 
lengths away from the atom (kgR 3> 1). This simplifies 
the theoretical treatments, and experimental approach 
greatly, yet allowing strong QED effects to be observed. 
Let us also mention here also that, since we use a single 
mirror, the QED effects for an exited atom are here easier 
to understand as a modulation of the electron coupling to 
certain modes rather than a change in the mode density 
0- 



I. ATOMIC DECAY AND LEVEL SHIFTS 




Region IL 



FIG. 1: Schematics of the mirror, and notations used to cal- 
culate the QED properties for an atomic dipole d in region I. 
The mirror has an aperture a and a radius of curvature R. 

It is well known that coupling an initially excited atom 
to a reservoir of electromagnetic field modes in the vac- 
uum state yields a spontaneous decay to the ground state 
and shifts its excited state energy. Also, when the atom is 
in the ground state, its energy is shifted due to absorption 
and re-emission of virtual photons, the so-called Lamb 
shift. When the electromagnetic field modes are modi- 
fied by nearby dielectric boundaries, these QED proper- 
ties are therefore also changed. Another picture, using 
radiation self-reaction only, can also be employed to de- 
scribe the modification of QED properties [35]. To find 
the relative contribution of both vacuum and self-reaction 
mechanisms, the dynamics of the corresponding quanti- 
ties in the differential equation one wishes to interpret, 
has to be Hermitian [23] 03]. When this is done, both 
effects are found to contribute to the same amount. 
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The free part of the light-atom Hamiltonian is the sum 
of the atomic and optical rest energies 



H = -hjJoa z 



E^ 



(1) 



where a z = — a 22 is the population difference between 
the upper and lower atomic states, and the creation 
operator for a photon in a mode \i of the reservoir, ujq is 
the atomic transition frequency, uj^ the frequency of the 
optical mode /i. The frequency of the optical modes is 
quantized using the boundary condition on a large cav- 
ity. In spherical coordinates we use a large quantization 
sphere of radius A, as depicted in Fig. [T] 

The interaction Hamiltonian in the Coulomb gauge 
and in the electric dipole approximation is 



Hmt = A{r,t) -p- 

mc 



2mc' 



:A 2 (f,t), 



(2) 



where p is the momentum of the atomic electron, r its 
position (i.e the position of the atomic nucleus in the 
electric dipole approximation) and m its mass, p will 
be written in terms of the electric dipole matrix element 
d of the two level atom as muiod/e x a y , where a y := 
i(o\2 — CT21). A is the potential vector. We decompose it 
over a complete mode basis e M as 



Am = E 



l2TThc 2 



e M (r) a^t) + h.c, (3) 



The excited state shift is due to the emission of real 
photons, whereas the ground state level shift arises from 
absorption and emission of virtual photons from the vac- 
uum reservoir [32]. To obtain the correct Lamb shift in 
the non-relativistic theory however, we have to add the 
term proportional to A 2 in the interaction Hamiltonian, 
which was actually discarded when getting to Eq. (JSJ. 
We come back to this in section [IV B| where we calculate 
the Lamb shift in the three-dimensional case. 

We would like to emphasize that as we only consider 
the atom to be far from the mirror (koR 3> 1, where 
fco = wo AO, the modification of the vacuum mode den- 
sity, affects the atom only if it is in the ground state [IS] . 
The change in the excited state properties are here due 
to a pure self- interference of the electromagnetic modes 
that the excited atom can couple to [32] [45] . We will 
then neglect the dependence of the coupling strength on 
Wfc for the excited level shift, which comes down to ig- 
noring energy level shifts (Van der Waals shifts) that are 
significant only if the mirror is very close to the atom. 
For a two level atom, these shifts are identical for the 
excited and ground state [IH] ; so the total 'near field' en- 
ergy shift would remain unaffected by the presence of the 
mirror anyway. 

Before moving to the three-dimensional results, and 
to gain physical insight on the process of emission and 
energy shifts in the presence of a mirror, we first perform 
a one-dimensional calculation. 



where the sum is to be taken over all normalized eigen- 
functions e M of the Helmholtz equation. 

In the Markov regime, that is when the reservoir and 
atom are correlated within a short time only, one can de- 
fine spontaneous emission rate and level shifts and, after 
solving the Hcisenberg equation, get to 



d_ 



(*(*)) - ~h(?)/2 + i (A e (r)-A g (r))](a(t)). 



(.) denotes the expectation value over a separable initial 
atom/vacuum-field state. We write 
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(r) = 2^| SM (rl 2 %- W o), 



1 



Un ±UJ Q 



(4) 
(5) 



where the +/— hold for the ground/excited state shift 
respectively, and 



9» (f) = W <M 



/ 2tt 

tuj,. 



d- e M (r)], 



(6) 



is the vacuum Rabi frequency of the mode 11. The sum 
over the /j, eigenmodes runs up to the Bethe momentum 
cut-off, K = mc/H [3J]. 



II. ONE-DIMENSIONAL MODEL 

In this section, we here assume that the atom can only 
couple to the electromagnetic fields through a set of one- 
dimensional spatial modes k along the mirror axis z This 
means that the other modes do not contribute to spon- 
taneous emission. This situation of course bears a simi- 
larity to the spherical mirror case, where also half of the 
light field is reflected back to the atom, and is therefore 
worth investigating in some details. 



A. Normal mode and quantization 

The scalar, one-dimensional, mode functions &k(z) of 
the problem must satisfy the Helmholtz equation 



V 2 e k (z) + \k\ 2 e k (z)=0. 



(7) 



The solutions are superpositions of plane waves trav- 
eling in reverse directions. We assume that the (one- 
dimensional) quantization domain has a length L, that 
the atom is at z — 0, and the mirror at z = —R. The 
density of k modes is readily found to be L/2-7T from the 
periodic boundary conditions at z — L. From the bound- 
ary condition on the mirror e k (z — —R) — 0, and after 
normalizing each mode to unity, we find that the mode 
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functions can be written 



V2LA 



(e lkz - e 



— 2ikR_ — ikz 



(8) 



where A is the transversal cross section area of the field 
around the atom. It is clear from this relation, that the 
wave going to the mirror (second term), and the wave 
going directly to free space (first term) will interfere per- 
fectly provided temporal coherence is fulfilled. 

Using this mode function, we now compute the influ- 
ence of the mirror on an excited atom, i.e on real photon 
processes . 



B. Real photon processes, atom in the excited state 

We here neglect the modification of the mode spectral 
density that couples to the excited atom, and consider the 
modification of the atomic state due to self-interference, 
as already discussed, so we neglect the dependence of 
on uj^. In the Markov regime, using Eq. Q and ([5]), we 
find that the atomic coherence decay ana excitedstate 
level shift at R = are then 

7(0) = 7 Fs(l-cos(2fc i?)) (9) 
A e (0) = 7FS sin(2fcoi?), (10) 

where the free space ID spontaneous emission rate is 



7FS 



2d 2 w 
HA ' 



(11) 



For 2kgR — 2im, (n an integer number), we get a com- 
plete suppression of spontaneous emission and excited 
level shift. On the other hand, for 2fc i? = irn, the spon- 
taneous emission is enhanced by a factor of two. 

The reason for such large effects is that the fields going 
to the mirror and the 'direct' fields can fully interfere in 
the Markov regime. To find out in which regime tempo- 
ral coherence is not-satisfied, so that the visibility of the 
single photon interference is not perfect, we would need 
to consider the exact quantum dynamical evolution with- 
out making a Markov approximation. Such an analysis 
would reveal that temporal coherence is reduced when 
the mirror is placed such that the time it takes for the 
light field to go to the mirror and back is larger than the 
atom decay time. This scenario was investigated, for ex- 
ample in [14] . and studied experimentally in [33], where 
signatures of non-Markovian dynamics were analyzed us- 
ing a Hanbury Brown and Twiss set-up. It was noted 
that bunching appears for short time scales, similar than 
the bunching that would appear for two classical sources 
(which here would be the atom and its far-distant mirror 
image) . 

In the 'extreme' non-Markovian regime, where the mir- 
ror is placed far away from the atom, a multimode field 
with width 1 /7fs is emitted towards the mirror with 50% 
probability. The other half goes to free space. The atom 
is completely in the ground state when this field returns 



from the mirror. It will re-excite the atom after a time 
t = 2R/c, but only partially since its temporal enve- 
lope is not the time-reversed spontaneously emitted field 
[To] and that its amplitude is twice as small. After such a 
(partial) re-excitation, another field will be emitted along 
the mirror so that the atom will again be re-exited later 
on. Eventually, the photon will leave after a complex 
dynamical process that resembles that of a multi-mode 
cavity [HI E2] . In this paper, we always assume a Marko- 
vian dynamics, where a linear decay and level-shifts can 
be defined. We will only be concerned with spatial deco- 
herence, by assuming ideal temporal overlap of the single 
photon with itself, i.e r <C 1/7fs- 

We assumed here that the atom only couples to the lon- 
gitudinal modes along the mirror axis, which is not real- 
istic for an atom in free space. We will now calculate the 
effects of polarization and use a spherical mirror. We will 
show that similar behavior appears for a full hemisphere, 
i.e in the Markovian regime, spontaneous emission and 
excited state level shifts can be suppressed. We also cal- 
culate the far field Casimir-Polder shift and compare it 
with the well-known calculations of the Casimir-Polder 
shift for an atom close to an infinite plane mirror. 



III. THREE-DIMENSIONAL MODEL 

In this section, the space will be divided into a part 
inside a sphere of radius R (region I), and the annular 
region between the sphere of radius R and the quantiza- 
tion sphere (region II). See Fig. [I] 



A. Normal modes 

In spherical coordinates, \i — (I, m, a) where a denotes 
the TE or TM modes, and (I, m) the quantum numbers 
for the angular momentum and spin respectively. The 
solution of the Maxwell equation for the electric field can 
be written as a superposition of the electric and magnetic 
multipoles [48, 



eTM(r) 



gi(k^r)Xi !m (ft) 

T^V X (/i(* M r)An>$))> 



(12) 
(13) 



where Xi tm (fl) = LY^ m ((l)/ + 1) is the vectorial 
spherical harmonic and // , gi are superpositions of spher- 
ical Bcsscl or Hankcl functions. Vt is the vectorial solid 
angle along the radial direction r. L = r x p/i is the 
angular momentum operator. The magnetic induction B 
is a superposition of the two multipoles 



(r) 



-V x ^(Lr^mtfl)). 



& TE (r) = fi(k^r)X ltm {n). 



(14) 
(15) 
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The radial functions are written as 

9i{k^r) = ciii{k^r) in region I, (16) 
= aihpik^r) + kj^k^r) in region II, (17) 

where ft, (fc^r) is the spherical Hankel function, and 
jl (kfj,r) the spherical Bessel function of the first kind. It is 
clear that, in the absence of the mirror, the 6; modes are 
the vacuum field from region II. We will show next that, 
in the single mirror geometry, we only need to quantize 
these modes. 



B. Quantization 

We will not solve the full quantum mechanical problem 
and quantize the electromagnetic field in the presence of 
the spherical mirror. This can be done exactly in the 
case of a full sphere [37] by solving the eigenvalue equa- 
tion derived from the boundary conditions. To the best 
of our knowledge, such a treatment has not been done 
for an opened geometry such as a hemispherical mirror. 
As was shown in [42] , the problem is not so difficult how- 
ever, if one assumes that the boundary condition lies in 
the far field (koR ^> 1), as we assume here, so that the 
mode density is close to that of free space. We are then 
mostly dealing with a continuum of modes, like in the ID 
calculations. 

Using the boundary condition on the large quantiza- 
tion sphere, we find that 



LA = l- 



(18) 



The density of free space vacuum modes is then 2A/n. 
From a point close to the focal point, the bi modes are 
non-degenerate if the mirror is large enough, so they are 
all orthogonal to each other. We normalize them so that 
each of them contains one photon. For the magnetic 
multipole we then require 



r 2 dr\ 9l {kr)\ 2 / dfi|X,, m (fi)| 2 



1. 

In free space, ai = 0, we thus get the condition 

A- 2 



\b,\ 2 



2A' 



(19) 



(20) 



where we used the fact that the main contribution to 
the vacuum fluctuations stems from the far field. The 
same relation holds for the vacuum modes of the electric 
multipole. 



region I, and associate the eigenmodes 6; to the contin- 
uous set k, so that — > A/n J dkJ2i m - F° r a dipole 
oriented along f for example, we can check that we get 
the usual free space spontaneous decay. Using Eq. Q, 
the formula for the field component along the radial di- 
rection (see Appendix), and setting a/ to zero in Eq.([l6|, 
we indeed find 



7fs - £ J dk^hm 



2-KUj 2 d 2 r . 
: — t-^ — d(u k - luq) 



n J; 2 (fe r ) ly |2 

1} w 1 l ' ml 



d 2 LO$ 

3-Khc 3 ' 



(21) 

(22) 
(23) 



which is the standard spontaneous decay rate of an atom 
in free space. In the last step, we used the addition for- 
mula for spherical harmonics J2 m l^,m| 2 = + l)/47r 
and the addition formula 



]TZ(Z + 1)(2/ + 1) 



1=1 



(kr)'- 



(24) 



One can show that the spontaneous decay rate is the 
same for a tangential dipole orientation. 

In the next sections, we will see how the mirror im- 
poses a fixed phase relation between the even and odd I 
modes appearing in the sum Eq. ( |24[ ) used in the sponta- 
neous decay calculation, and how this modifies it. This 
is already hinted by noting that 



£ 

/ even/odd 
rr/2 



Z(Z + 1)(2Z + 1) 



(kr) 2 



d6 sin J 9 



sin 2 (fcr cos 9) 
cos 2 (kr cos 9) 



(25) 



where the odd/even modes correspond to the sine/cosine 
functions. Depending on whether the atom is at the node 
or the antinode of the standing wave, it couples preferen- 
tially to the even or odd modes. Without a defined phase 
relation between even and odd / modes, as is the case for 
free space, their spectral densities always adds up to 2/3 
as we just saw, but they can otherwise cancel or add up 
coherently. This will allow suppression or enhancement 
of the even or odd vacuum modes fluctuations, and thus 
significant changes in the decay and shifts. 



D. Boundary conditions 



C. Free-space decay 

Having normalised the vacuum modes bi , and the nor- 
mal modes of the system, we can calculate the distri- 
bution of vacuum fluctuations and atomic properties in 



The boundary conditions on the electric and magnetic 
fields are 



n x E l \ r=R = nx E u \ r =R, 



ft x B l 



r=R 



ft x B u \ r=R , 



(26) 
(27) 



() 



for 9 — [a,ir[ and (j> = [0, 2tt[ (see Fig. 1), where n is 
the normal to the mirror. We write £r ,n , the electric 
field in region I, II. Assuming that the mirror is a perfect 
conductor, as we always do in the paper, we also have 



n x E l \ r = R 



n-B l 



r=R 



= o, 
= o, 



(28) 
(29) 



for 8 = [0,a[ and 4> = [0,2tt]. 

Two sets of equations for the transverse electric and 
magnetic multipoles can then be obtained and solved for 
ci to calculate the total field in region I. 



E. System of equations 

We use the symmetry along <p by projecting the bound- 
ary conditions over the m modes (multiplication by 
and integration over 4>). From equation ( [28] ), we get, for 
example for the 8 component of TM mode, the relation 

oo 

J29i(kR)Ai, m (6) = 0; for 6 = [0,a[, (30) 



/=i 



where 



V (t + m)! + 1) 

and P™ is the associated Legendre polynomial. After 
projecting over m and I, and using the orthogonality of 
the spherical modes, we require 

ciii{kR) = ai h ( p(kR) + hji{kR). (32) 
Last, we use Eq. |27|) and obtain 



Or 



for 9 = [a, 7r], and r = R. Using the Wronskian for 
spherical Bessel functions, we then get the two sets of 
equations 

oo . 

E( c <-M^t^ = °; * = M 04) 

h\ '(kR) 

OO 

J2 c tMkR)Ai,m = 0; 0=[a,7r] (35) 



i=i 



i=i 



for the coefficients of the magnetic multipole. Similarly, 
we get 

A 



£ (d '~ e;) r J)Z , v = °i = [°'«[ ( 36 ) 
z=i [rftj '(At)]' 

OO 

]>>[nKfcr)MU = 0; 6=[a,ir] (37) 



for the electric multipole, where [rji(kr)]' = 
d{rji(kr))/dr\r = R, and A' lm = dA Lm /d9. e, 
and di are the amplitude coefficients for fi(kr) equiva- 
lent to ci and 6; for gi(kr). Each set of equations can 
be written as a Fredholm equation that can be solved 
numerically [56] . 

Here, we will show that an analytical solution can be 
found in the limit where k$R 3> 1(1 + 1), that is using a 
large mirror and/or looking at the field fluctuations close 
to the focus. We note that the field is orthogonal to n 



far from the origin (See Eq. (68)). The condition (29) is 
then satisfied automatically. In this case, the solutions 
of the two sets of equations will be identical for both the 
TE and TM modes, so that the reflection off the mirror 
will preserve polarisation. 

The solution can then be found using scalar fields. We 
will include the polarisation dependence of the dipole 
emission later, after having identified far field plane wave 
modes. We will here write 



b b (f) = Y*jt(kr)Yt 



l,m-i 



(38) 



1.771 



for the total field amplitude in region I, and solve for q 
as a function of b\ using the two equations 

oo 

= 0; 6=[0,a[ (39) 



1=0 



h ( , L >(kR) 



J2cijl(kR)Yi, m = 0; 9=[a,n}. (40) 



1=0 



We simply removed the 1(1 + 1) dependance of the field 
modes, which as can be found from Eq. (25), is equiva- 



lent to ignoring the polarization dependance of the dipole 
emission (sin 2 9 in the integral) . Note that the summa- 
tion starts at I — now, as is allowed for scalar fields. 
This is not true for the general solution of the (vecto- 
rial) Maxwell equation which does not have spherically 
symmetric solutions. 



F. Full-hemispherical mirror 

We first assume that the mirror covers 2tt steradian, so 
a = it/2. Furthermore, since kgR 3> 1(1 + 1) we expand 
the Bessel functions in the far field. 



^ (1) (fci?) w (-i) l+L e lkH /kR, 



(41) 



ji(kR) w sm(kR-lTr/2)/kR. (42) 
Using the relation 



i=i 



(43) 
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in Eq. ( 39 ) and Eq. ( 40 1 and then solving for q , we obtain 



after some algebra, the two relations 



ci> 



JkR 



COS 

Il',l,m) for /' even, 



ikR , 



cv = ie"'"-' sin(fci?)(6;/ + i £ for I' odd(44) 

where 



I'l.m 



(-1) 



(Z+i' + l)/2 



<(, J6=[0,ir/2] 



Y Wm Y^ m dn. (45) 



We denoted J2 a e ^ ne sum over odd/even I modes. The 
two terms on the right hand side of each equation are the 
reflected and incoming vacuum amplitudes contributions 
to the field in region I. 

At the focus, only the V = mode is predominant since 
the radial amplitudes (given by jv(kr)) are negligible for 
higher The total field <t>b{r) at the center will then be 
zero for kR = nir, as can be seen from Eq. (44) and (38 1. 



If we draw any standing wave from the mirror through 
the origin to the other free space boundary with the con- 
dition that kR — mr, it is then invariant under rotation 
upon 9. It follows that to completely describe the field 
in region I when the total field is zero at the focus, we 
only need even modes (as is apparent from Eq. (44 1). 
For the case where kR = nir/ 2, we will however need the 
odd modes to find the total field in region I. This two 
extremal conditions explain why even and odd I' modes 
behave differently with respect to mirror positions, as we 
already anticipated in section ( III C I . 



We also note that, contrary to what one would find for 
a mirror covering 47r (which would behave as a cavity), 
the angular asymmetry of the hemisphere does not yield a 
one to one mapping of the free space modes b[ of region II, 
to the c\ modes of region I. The total far field amplitude 
that can enter region I is here necessarily a superposition 
of even and odd modes. 

We already showed that the total field at the focal 
point <pb (r = 0), will be zero at a node. Spontaneous 
emission and level shifts will then certainly cancel. To 
demonstrate this result, and also treat the case of a fi- 
nite mirror size, it is useful to introduce some relations 
between the plane waves and spherical harmonics. 



G. Plane wave decomposition 



Let us first write the far field amplitude in region I as 



At 



where 



m,e/o 



(47) 



is the scattering amplitude for the even/odd mode. It is 
easy to show that the superpositions 



rm = - 2 {f°-n 



n\n) = -(f° + n, (48) 



correspond to incoming and outgoing angular amplitudes 
respectively, and that they are connected via 



(49) 



Here P is the parity operator. This relation shows the 
Gouy phase shift acquired by the incoming field as it 
turns after focussing into an outgoing field. The same 
relations also hold for the far field amplitude g(Q) of the 
field in region I, where one can define 



m,e/o 



(50) 



The total field amplitude at any point of region I can in 
fact be written as a superposition of plane waves weighted 
by the far field amplitudes g(Ct). Using the expansion of 
the plane waves in spherical harmonics 



£ i\2l + l) ]l {kr)Y u 

i,(e/o) 



cos(fc ■ r)' 
sin(fc • r) 



(51) 



and Eqs. (50) and (38), we indeed obtain 



i 

2^ 



2n 



d(l[g e cos(fc • r) + ig° sin(fc • r)} ,(52) 



where we wrote := ^ f e=[Q /2] for simplicity. The 
field inside region I, is uniquely given by the coefficients 
g e /° : which are linked to the incoming vacuum modes 
/ m via the boundary conditions. This treatment was 
also used in [42] for an open cavity. 



H. 



Finite size mirror 



To treat the case of finite mirror size, we decompose 
the boundary conditions into three interfaces: the mirror 
(9 = [0, a[), the opposite free space interface between 
region I and II (9 — [tt — a,7r]), and the remaining free 
spac e (9 = \a, n—a[). As before, we project the equations 
Eq. d39]40l) over the m modes and using the results of 
section ( IIIG| ) get to the results 



(l - e 2lkR T(9))g° = -2e lkR sm{kR) (f in + T{9)f out ) , 
(l + e 2ikR T{6))g e = -2ie lkR cos(kR)(f in - T(6»)/ out ), 

(53) 

where T{9) is a function that is zero for 9 = [0,a[, and 
unity for 9 = [a,7r/2[. One can check that setting T{9) 
to zero for 9 = [0,7r/2], and going back to the spherical 
basis, we recover Eq. (44). 



Using Eq. (49), and the fact that T is only zero or one, 



we can rewrite Eq. ( 53 1 as 



T{6)(1 -P)+ 2i{\ - T{9))e lkH sin(kR) f 



T(0)(1 + P)- 2(1 - T(9))e lkH cos(fei?) 



r, 

(54) 



which uniquely relates the incoming field from region II 
to the far field modes in region I. With the use of the 
relation ( 52 ) , we finally get to the result 



[.ill 



(55) 



where 



T, 



2 1 



dn r 

4?r 



T(6){e lk - r + e- tk - r P) 



+ 2(1 - T(6))e lkH cos(fci? + k ■ r ) 



Y, 



(56) 



The total field in I is a superposition of waves that 
travel through the point of coordinate r without being 
reflected by the mirror (first term in Eq. 56 1 and waves 
that are reflected by the mirror (second term) as one 
could have easily figured out. A similar result was found 
for an open cavity with variable reflectivity in |42j , where 
the density of modes of a single mirror was calculated 
by setting one of the mirror reflectivities to zero. Here 
we derived the normal modes amplitudes for the case of 
a single mirror by using a Neumann condition between 
regions I and II (through continuity of the B field) . 

This mode function can now be used to calculate the 
density of fluctuations around the focal point, to ob- 
tain the change in excited level shifts and spontaneous 
emission rates, and find the ground state Casimir-Poldcr 
shifts as a function of the mirror's numerical aperture. 



Atom at the antinode 



ro 1.4 



OJ 



1 

0.8 - 
0.6 - 
0.4 - 
0.2 

°o- 



Atom at the node 



.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Numerical aperture 

FIG. 2: Spontaneous decay rate as a function of the spherical 
mirror numerical aperture, for an atom at the mirror focus. 
Top curve : The atom is at the anti-node of the standing wave 
formed by its retro-reflected field. Bottom curve : Atom at 
the node of the standing wave. Large changes in the spon- 
taneous emission rate are already expected for moderate nu- 
merical apertures. 



The dependence of the function cf) on the b modes is now 
explicitly given by k. In the last equation, we used the 
fact that the spherical harmonics form an orthonormal 
set of modes. We also introduced p = 1 — T, the reflec- 
tivity of the mirror. 

We can now also include polarisation, and get 



l*e(*0l s 



k 2 


f d&3 


1 - 


d- (i 


2 i 


A 




d 





1 - pcos(2(kR + k ■ r)) 



(60) 



If we set p to zero everywhere, we recover the density of 
vacuum fluctuations in free space. 



IV. QED EFFECTS CLOSE TO THE FOCUS OF 
A SPHERICAL MIRROR 



We can now calculate the QED effects on the atomic 
electron from Eqs. Q, ^ and (56). The modified decay 
and level shifts are given by 



and 



, . 2d 2 Au 2 f dk, , _., n 



he Ik 



1 



k ± fcn 



(57) 



(58) 



where 



(1 - p(ft) cos(2(fci? + k-r)). (59) 



A. Real photon processes 

The spontaneous decay normalized to the free space 
decay rate 7fs, is 



7(f) 



d-Ct 



3 f dQ 
2 J 2. 4n 
1 - pcos(2(k R + k ■ f)) 



(61) 



This quantity is plotted in figure [2j as a function of the 
numerical aperture (defined as NA= sin(a)), for a lin- 
early polarized dipole positioned at r = 0, and orthogo- 
nal to the mirror axis. Spontaneous emission is found to 
vanish for a mirror position such that cos(2fcoi?) = 0, i.e 
at the node. A twofold increase in the spontaneous emis- 
sion rate is found when cos(2fc i?) = 1, at the antinode. 

We note that for the numerical aperture used in |12j . 
(NA=0.4, i.e 4% of solid angle), a spontaneous emission 
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rate change of 24% is predicted. Such a, perhaps un- 
expectedly, large change of the decay rate may be un- 
derstood by noting that the factor of two coming from 
the interference between the two reflected and direct am- 
plitudes, translates into a factor of 4 in density (for 
'small' numerical apertures). Together with the inclu- 
sion of the polarisation properties of the dipole emis- 
sion, another factor of 3/2 is gained, which in total gives 
(3/2 x 4) x 0.4 = 0.24. The difference between the ob- 
served 1% change of the excited state population in [12] 
and the 24% modification of the spontaneous emission 
predicted here can be partly explained by residual atomic 
motion, finite spatial overlap, finite temporal coherence 
or multi-level effects in the experiment. 

We can compute the modification of the decay rate as 
a function of distance from the focus to estimate the sen- 
sitivity to mirror, or lens, positioning. Fig. [3] shows the 
dependance of 7 for an atom that is displaced away from 
the focus of a hemispherical mirror and where the mirror 
is positioned such that |<^e(r = 0)| 2 = 0. The sponta- 
neous emission rate is close to zero within a volume A 3 
around the focus, and oscillates for a few wavelengths 
until it reaches the free space value. More precise formu- 
lae must however be used when the atom is far from the 
focus, as the approximation kR 1(1 + 1) is no longer 
valid for large distances from the focus [4*2"] . 



<D 
T3 




2 4 6 8 10 

Distance from the focus (in units of X) 

FIG. 3: Spontaneous decay rate of an atom as it is displaced 
from the focus, in the case of a full hemispherical mirror. 
Here the mirror is positioned such that there is a node of the 
standing wave at the focus. Trace (i) corresponds a scenario 
where the atom is displaced along mirror axis, (ii) is when 
the atom displaced perpendicularly to the mirror axis. 



The excited state level shift will also be altered in the 
same way due to a modified coupling to the retroreflected 
modes. We find after contour integration, that the ex- 
cited state shift, normalized to the free space decay rate, 



A e (r) 



2n r 

psm(2(k a R + k -r)) 



2 ./ 2w 4^ 



This gives an oscillatory level shift of amplitude p at 
the focus. For a full half-mirror, the level shift com- 
pletely cancels for mirror positions such that k$R = inn, 
where the decay rate is also zero. Its evolution with 
numerical aperture is similar to the spontaneous emis- 
sion rate change. Such large level shifts variations can 
yield a strong confining potential and would be interest- 
ing to study experimentally with large numerical aper- 
tures, similarly to what was done in (50] . 



B. Virtual photon processes 

The Lamb shift of the ground state can be computed 
in the very same way as the excited state shift. In the 
simple case of a full half-mirror, and with the atom at 
the mirror focus, we get 



A 9 (0) 



7fs 
£; 



c/fi 



dk 



fc + k 



sin 2 (kR). 



(63) 



We can write this result as a sum of three terms that 
can be easily integrated : The electron self-energy, the 
free space Lamb shift, and the Casimir-Polder shift. The 
electron self-energy is 



K = 7FS 



K 
2ko 



1 



sin(2KR) 
2KR 



(64) 



where we wrote K — mc/H. This quantity is identified 
by setting ko to zero in Eq. ( 63 1 and in fact cancels with 



the shift from the A 2 part of the Hamiltonian, as can be 
easily checked. This procedure is known as mass renor- 
malisation [32]. The free space Lamb shift is 



A{ s = 7fs log 



(65) 



The modified Lamb shift (or Casimir-Polder shift), the 
only observable quantity, is 



= 7FS 



2k ° R dx 



2KR 



cos(x — 2k R). 



(66) 



The Casimir-Polder shift goes to zero for very large 
mirror-atom distances (2k^R 3> 1) as expected, and can 
be approximated by 



7FS 

(k R) 2 



(67) 



closer to the mirror (but always at the focus). We note 
that Ag P drops slower with distance than in the plane 
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mirror case, where it decreases like 7fs/(^o-R) 4 [27]. The 
difference lies in the fact that there, the mirror does not 
cancel as many electromagnetic modes, which yields a 
faster decay of the QED effects with distance. 

As an example, using a decay rate of 15 MHz, a wave- 
length of 493nm (the Si/ 2 to Px/2 transition of Barium 
138Ba + ) and a mirror distance of 1 cm, gives a Lamb 
shift of 100 Hz which is experimentally measurable us- 
ing modern spectroscopic tools [53]. The complete level 
structure will have to be used for a precise estimation 
of the total shift [3U 0S] , but this result is encouraging 
for experimental investigations of Casimir-Polder shifts 
using trapped atoms. 

The associated force on a trapped atom due the 
Casimir-Polder shift may be computed easily from 
the above formula [35]. It would be interesting to 
study how Casimir photons created by an oscillat- 
ing mirror |55j . and also real photons 'modulated' 
by the mirror |54j . affect the center of mass motion 
of a trapped atom with a high numerical aperture mirror. 

To summarize these results, the modified decay and 
both ground and excited level shifts are plotted in Fig. [4] 
for a full hemi-spherical mirror. The decay rate shows 
undamped oscillations between full suppression and max- 
imum enhancement as a function of mirror distance (or 
mirror radius of curvature). The excited state level shift 
oscillates between —7 and 7, whereas the ground state 
shift damps out as a function of mirror distance. 



-1 













Lamb shift 


Excited state 



5 10 15 20 25 30 

Mirror distance from the atom (in units of X) 



FIG. 4: Normalized spontaneous decay rate and excited and 
ground state shifts of an atom as a function of the mirror- 
atom distance, in the case of a full hemispherical mirror. 



suppress the atom's spontaneous emission and excited 
level shift. We first presented a one-dimensional treat- 
ment that explained the underlying physics behind the 
full spherical mirror scenario. The modification of QED 
atomic properties was then calculated as a function of the 
spherical mirror's numerical aperture beyond the parax- 
ial approximation. Large effects are found for moderate 
numerical apertures, and with mirror-atom distances of 
several wavelengths when the atom is located at the mir- 
ror focus. 

This result is also relevant for the growing field of free 
space coupling to single absorbers, where full absorption 
of a single photon field requires a large coverage of the 
spatial dipole emission with the incoming spatial mode. 
The single hemi-spherical mirror system may here serve 
as an efficient quantum memory that can release a stored 
excitation on demand on a two-level atom transition, by 
controlling the mirror position in a dynamical fashion. 
As an application of our calculations, one also expects 
full super/sub radiance with two atoms interacting via 
large lenses covering only half of their respective dipole 
emission profiles, see for example [23 [52] for studies of 
this effect. Last, atom trapping using the dipole force can 
be very efficient here, due to the steep spatial dependence 
of the level shift across the atom. 

Finally, we calculated the Lamb shift, and showed a fa- 
vorable scaling of the spherical geometry over the plane 
mirror case. The Lamb shift scales as 7fs/ (koR) 2 , where 
R is the mirror radius of curvature. This contrasts with 
the 7Fs/(fco-R) 4 scaling law found for a plane mirror. 
Using Rydberg atoms and high numerical aperture el- 
ements can potentially yield very large shifts even for 
atom-mirror distances of a few centimeters and serve as 
a precise test-bed for investigations of Quantum Electro- 
Dynamics. 
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Appendix Along 9, we have 



We can decompose the field in region I into radial and 
longitudinal parts using the relation 



tm(0 = ZZT7>9i{r)— ^!; m = (71) 



k kr or \ r> F) Y, 

+ (68) - (72) 



where ft = r/\r\. Using the decomposition (68), and 

the eigenvalues of the angular momentum operator L, And along r 

we then get an expression of the electric and magnetic 

field in spherical coordinates. 

The electric field multipoles along read 



a Y , e T TE (f) = -i fi{r)Y hm . (73) 

U 1 l.m Sr>r\\ rZT 



e „,(n = -mlr) ¥e7 = (69) 

= ra^Jfrj <™» 
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